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Abstract 
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^ ] New exact solutions of Einstein's gravity coupled to a self-interacting 

Pfj' conformal scalar field are derived in this work. Our approach extends 

a solution-generating technique originally introduced by Bekenstein for 
massless conformal scalar fields. Solutions are obtained for a Friedmann- 



$H ' Robertson- Walker geometry both for the cases of zero and non-zero cur- 

vatures, and a variety of interesting features are found. It is shown that 
one class of solutions tends asymptotically to a power-law inflationary be- 
haviour S{t) ~ t^ with p > 1, while another class exhibits a late time 
approach to the S{t) ~ t behaviour of the coasting models. Bouncing 
models which avoid an initial singularity are also obtained. A general dis- 
cussion of the asymptotic behaviour and of the possibility of occurrence of 
inflation is provided. 
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1 Introduction 

The action describing the interaction between the gravitational field and a 
scalar field is generally constructed using the minimal coupling principle. Only 
such terms are included that give the curved generalization of the fiat Minkowski 
spacetime form of the equation of motion for the scalar field, i.e., of the Klein- 
Gordon equation. Terms that contain the spacetime curvature are usually disre- 
garded. 

Recently, a non-minimal coupling between the spacetime curvature and a 
scalar field has attracted a great deal of attention as a possible improvement of 
the minimal coupling case [1-10]. In principle, the consideration of such an in- 
teraction allows one to take more properly into account the infiuence that, in the 
very early Universe, the extremely high value of the curvature potentially had 
in the dynamical behaviour of this coupled system. This point of view seems 
reinforced when the low energy limit of unified theories (e.g., Superstrings, Su- 
pergravity, Kaluza-Klein theories,. . . ) is considered. Upon compactification to 
four spacetime dimensions, the effective action obtained in those theories often 
exhibits the aforementioned feature [11]. 

If one requires that such a non-minimal coupling should: (i) involve only the 
scalar field and not its derivatives, and (ii) be characterised by a dimensionless 
coupling parameter; one is led uniquely to the following action [1,2,12] 
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where g is the determinant of the metric gab, Mpi = Gj^ is the Planck mass, 
R is the Ricci scalar curvature, is a real scalar field singlet, V{(j)) is a general 
(effective) potential expressing the self-interaction of and ,^ is a dimensionless 
parameter that characterises the strength of this finite coupling between the scalar 
field and the curvature scalar. Note that the 1^0^ i? term in the action (1) is the 
only conceivable local scalar coupling satisfying the requirements (i), (ii) stated 
above [12]. (Our sign conventions are that gab has signature (— , -|-, -|-, -|-) and the 
Riemann and Ricci tensors are defined as i?^^^ = ScF^^ — • • • , Rab = Rad,- ^^ ^^^ 



units such that h = c = 1, and usually we shall set 8n /Mpi = SttGn = 1 for the 
sake of convenience.) 

For ,^ = we recover from eq.(l) the usual action for a minimally coupled self- 
interacting scalar field. On the other hand, for ,^ = | we obtain the conformal 
coupling case, and if V{(j)) does not contain dimensional scales (i.e., the potential 
is zero or a quartic power of 0), it can be shown that the Klein-Gordon equation 
is invariant under a conformal rescaling of the metric and the field itself, i.e., is 
conformally invariant [13]. 

From expression (1), we find that its possible to define an effective Newton 
gravitational constant G^ , which is given by the inverse of the coefficient of the 
curvature scalar, depending on the scalar field as 

^'^^ = M2,-87r02- (2) 

Thus, if the scalar field is assumed homogeneous, the essential role of the non- 
minimal coupling is to induce a time variation of the Newton gravitational con- 
stant. 

The theory based on (1) has been studied by a number of authors in several 
different contexts. Hosotani [1], has discussed how the back-reaction of gravity 
affects the stability of a scalar field, thus clarifying the conditions for this theory 
to have an absolutely stable ground state in which is constant and the spacetime 
background is maximally symmetric (in connection with this issue see also Refs. 
[2] and [3]). Dolgov and Ford [4], in turn, have investigated the possibility of a 
damping of the cosmological constant, in an attempt to solve the cosmological 
constant problem. In Ref. [5] the canonical quantization procedure is applied to 
(1) and a homogeneous and isotropic minisuperspace model is discussed. Finally, 
the prospects and consequences of infiation within this theory have been consid- 
ered in a variety of works [6-8]. However, there aren't many exact solutions in 
this context [9,10,14]. 

Here, we aim at addressing this question in the framework of the confor- 
mal coupling case. We obtain new exact solutions of the coupled Einstein-self- 
interacting conformal scalar field equations in the particular case of a spatially 



homogeneous field propagating both in fiat and curved Friedmann-Robertson- 
Walker (FRW) spacetime bacligrounds. This is of considerable interest to assess 
the cosmological implications of this theory and is done by extending a solution- 
generating technique originally introduced by Bekenstein for massless conformal 
scalar fields, which is based on the use of a conformal transformation [14]. This 
device, discussed at length in Refs. [8,14,15], enables one to obtain a representa- 
tion of the theory in which the original coupled system is described by Einstein's 
gravity plus a redefined self-interacting scalar field, which is now a minimally 
coupled one. This new frame, sometimes called the Einstein frame, is particu- 
larly suitable to the search of new exact solutions (e.g., following the strategy 
recently adopted by Ellis and Madsen [16]). Given a solution in this new frame, 
the extended solution-generating technique allows one to get a solution in the 
physical frame. 

This paper is organized as follows. In the next Section, we derive and review 
the general field equations for the non-minimal coupling model characterised by 
the action (1). We also present here the equations that govern the evolution of 
a FRW spacetime in the presence of a homogeneous self-interacting conformal 
scalar field. In Section 3, we present the extended solution-generating technique, 
reviewing in some detail the conformal transformation upon which it is build and 
emphasizing its usefulness regarding the derivation of exact conformal scalar field 
cosmologies. In Section 4, we apply the extended solution-generating technique to 
a number of classes of solutions in the Einstein frame that were recently obtained 
by Ellis and Madsen [16]. The classes of solutions obtained in the physical frame 
exhibit a diverse range of properties. For instance, one class of solutions tends 
asymptotically to a power-law inflationary behaviour, while another class exhibits 
a late time approach to the S{t) ~ t behaviour of the coasting models. Bouncing 
models which avoid an initial singularity are also presented. Finally, we conclude 
with a Section devoted to an overall summary and discussion of the results. 

2 The Field Equations 



We begin this section by deriving and reviewing the general field equations 
for the non-minimal coupling model characterised by the action (1). 

Variation of the action (1) with respect to the gravitational degrees of freedom 
yields Einstein's equations 

Rab — -QabR = Tab, (3) 

where Rab and Tab are the usual Ricci and energy-momentum tensors. The later 
has the form 

Tab = {l-C<PY'[{l-20da<pdb<l)+{2^-^)gabdc<pd'<P-2^(j)Va{db(P)- 

gabVi(P) + 2^gab(pa<pi (4) 

with Va standing for the covariant spacetime derivative associated with the met- 
ric Qab and □ = VqV" is the corresponding curved spacetime D'Alembertian. 
Variation with respect to the scalar field results in the Klein-Gordon equation 

□0_^0i?_V'(0) = O, (5) 

where V'{(f)) = dV{(j))/d(l). 

Next, we specialize the field equations for the particular case of a FRW space- 
time background filled with a homogeneous self-interacting conformal scalar field. 

The line element for a FRW spacetime can be written in standard normalized 
comoving coordinates (t, r, 6, ip) as 

ds^ = -df + S\t) i^K-^ + rHde^ + sin^ Od^^)] , (6) 

yl — kr^ I 

where Sii) represents the cosmic scale factor and k is the curvature index [k, = 
0, ±1 corresponds to the fiat, closed and open FRW model, respectively). The 
independent components of Einstein's equations (3)-(4) for this metric can be 
represented by the Friedmann and Raychaudhuri equations. In the conformal 
coupling case they become 

3H^ + 3K = {1- CV')"' [3C'0' + V^(0) + 00^1 , (7) 

3H + 3H^ = (1 - C(j)'^y^ [V(0) - 3C^(0^ - 00if - 



where, in accordance with the FRW symmetry, we have assumed that = 0(t). 
Here K is the purely spatial part of the scalar curvature K = k/S'^X = (1)^^^ 
and, as usual, H = S/S denotes the Hubble parameter and the overdot means 
derivative with respect to the synchronous cosmic time, t. On the other hand, 
using eqs.(7) and (8), the Klein-Gordon equation (5) becomes 

4> + 3H^ + C'0(1 - CV')"' [4V^(0) + 300if + 00] + V^'(0) = 0. (9) 

It is worthwhile to point here that, for the case of a homogeneous scalar field 
propagating in a spatially homogeneous and isotropic spacetime background, the 
energy-momentum tensor (4) can be rewritten in the form of a perfect fluid for 
any value of ^ [2,7]. The energy density and pressure of the equivalent perfect 
fluid are, in the conformal coupling case, given by 

p^ = {i- e<P'r'[K'<P' + v{<p) + HHi (10) 

and 

p^ = {l- CV)~'[C'(0' - 200 - 400/7) - F(0)], (11) 

respectively. It is evident that neither the energy density nor the pressure are 
sign-definite. As a consequence, \{p^ + 3^,^) has indeterminate sign, so that the 
strong energy condition [17] 

(T„, - \Tg,,)U'^U' = i(p^ + 3p^) > 0, (12) 



where f/" = d"'(f)/J—db(l)d''(j) denotes the 4- velocity vector of the equivalent perfect 
fluid, may or may not be valid. In the later case, the Hawking-Penrose theorem on 
the existence of singularities [18] is not applicable and singularity- free solutions 
to equations (7)-(9) are then possible. In Section 4 a number of singularity-free 
solutions of those equations will be presented. 

We will find it convenient to define a new scalar field variable -0 such that 
■0 = ((f). Using this variable eqs.(7)-(9) become respectively 

3i/2 + 3K = {1- V^')-'[3^2 + V{^) + C^ijijH], (13) 

3H + 3H^ = (1 - ipY'^lViip) - 3{iij^ - i^i^H - ^V^)], (14) 



i) + ?>Hi) + ^(1 - ^2)-i [4C V(7/') + 3^^if + ^^] + C^\/'(^) = 0, (15) 

where now V'{il)) = dV (ip) / dip . 

3 The Extended Solution-Generating Tech- 
nique 

As we remarked in the Introduction, the solution-generating technique em- 
ployed in this work is based on the use of a conformal transformation of the 
spacetime metric, which reduces the theory under consideration to one contain- 
ing canonical Einstein's gravity plus a canonical scalar field. In what follows we 
present this solution-generating technique, reviewing in some detail the conformal 
transformation upon which it is build. 

Let us perform the following Weyl conformal transformation 

gab = ^'9ai, , ~9''' = n-'9'' , ~g = n'g, (16) 

where 

n^ = \i-^<p% (17) 

Since the transformed scalar curvature is 

R = n-^[R-Qn-^nn], (is) 

it is then straightforward to show that action (1) can be written in the Einstein 
frame as 

S[r\(l>] = Jd'x^ 

where 

F\<j>) = 

and 

In order to put the kinetic term of the scalar field in the canonical form, we define 
a new field as 

0=/rf0F(0). (22) 



l- v(0)9„09>-\/(0) 


(19) 


1-^(1-600' 

(i-e0T ' 


(20) 


\/(0) 


(9.1) 



The theory can thus be written in terms of these new variables as a minimally 
coupled theory with the action 



I - i9,09> - V(0) 



(23) 



S[r',4>]=J d^x^-~g 

Taking ,^ = ^ in eq.(22) and solving for yields 

0(0)=rHanh(C0), (24) 

and upon substitution of this result in eq.(17) one obtains 

n-\i) = cosh(C0(t)), (25) 

and eq.(21) can be rewritten as 

n0) = (l-CV)V[^ln(i±|^)]. (26) 

From the point of view of the search of exact solutions to the field equations (3)- 
(5) in the conformal coupling case, this changing frame process can be interpreted 
in the following way: if the set {(jab, 0, V{(f))) forms a solution to the field equations 
in the Einstein frame, that is, a solution to eqs.(3)-(5) but with .^ = 0, then the set 
{gab, 0, V{(p)), where gab, <P and V{(f)) are given, respectively, by eqs.(16), (24) and 
(26), with Q given by eq.(25), constitutes a solution to the same field equations 
but in the physical frame. The massless case of this solution-generating technique 
was given by Bekenstein [14]. 

Note that this solution-generating technique is valid for any metric gab- In the 
particular case that concerns us here, gab is a FRW metric, whose line element 
was given in eq.(6). The line element of the conformal spacetime is, according to 
eq.(16), given by 

rf§2 = Q^ds^. (27) 

This can be brought into the same form as eq.(6) by the following change of 
variables 

di = Qdt, (28) 

S = QS. (29) 
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Then, we have in the Einstein frame 
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rf§2 ^ _^p ^ ^2(^~) ( ^^ + ^2.^^2 ^ gi^2 ^^^2^ ] ^ .ggN 

yl — /cr^ j 

i.e., the metric in the Einstein frame cjab is also of the FRW type, but with a 
redefined scale factor. 

Before proceeding further, it will be useful to pause in order to introduce some 
observational quantities of interest. The cosmological density parameter Vld (we 
adopt this notation slightly at variance with the standard usage to avoid any 
confusion with the conformal factor ^2 ) is defined by 

«.^^. (31) 

and when combined with the Friedmann equation (7), can be rewritten as 

n, = l + t. (32) 

The deceleration parameter q is defined by 

Following [19], we take the sign of the deceleration parameter as indicating 
whether a given cosmological model inflates or not. The positive sign corre- 
sponds to "standard" decelerating cosmological solutions whereas the negative 
sign indicates inflation. Although there are several definitions of what is actually 
meant by inflation, this is the most commonly used in the study of exact scalar 
field cosmologies. 

Defining a "mixed" cosmic scale factor by S{t) = S(t{t)) we have 

S = QS' (34) 

S = Q{Q'S' + QS"), (35) 

where the prime denotes differentiation with respect to the cosmic time t of the 
Einstein frame. A very convenient form for Q^ and q can then be obtained 
inserting these results into eqs.(32) and (33) 

fi. = 1 + ^, (36) 



__n's_ss" 

For future reference, we end this section giving the equations of the extended 
solution-generating technique for a spatially homogeneous and isotropic metric 
(eqs.(24)-(26),(28)-(29)), rewritten in terms of the ■?/' scalar field variable 

S{t) = cosh{i! {i))S{i), (38) 

^ = cosh(^(t)), (39) 

at 

i){t) = tanh(V^(t)), (40) 



vw = ii-ryv 



r\^/i+^^ 



(41) 



2 Vl-^. 
This equations can be looked upon as a map between the space of Einstein frame 

solutions and the space of physical frame solutions. 

4 The Solutions 

We now apply the solution-generating technique developed in the preceding 
section to a number of interesting situations. We consider, in the Einstein frame, 
the solutions which were recently obtained by Ellis and Madsen [16]. Here only 
three of the five exact solutions found by Ellis and Madsen will be mapped to the 
physical frame. 

The titles adopted in the following subsections refer to the behaviour of the 
solution in the Einstein frame. 

4.1 Power-Law Expansion 

First, let us consider the solution in the Einstein frame characterised by a 
power-law behaviour 

S{i) = Af\ (42) 

i!{i) = lniP, (43) 



y(V^) = \/oexpU-V^j, (44) 

where A,n are strictly positive constants with n 7^ 1, p{n) = (n/3)"'^/^, Vo(n) = 
n{3n—l), which is only possible for k = 0. This class of solutions is inflationary for 
n > 1, because if the power-law expansion continues for long enough, it will solve 
all the well known problems of the standard hot big-bang model (horizon, flatness, 
spectrum of primordial density perturbations,. . . )[20]. However, it behaves like 
a "standard" decelerating big-bang when < n < 1. 

Mapping the class of solutions (42)-(44) in the Einstein frame into the physical 
frame with the help of eqs.(38)-(41), we obtain the following parametric class of 
solutions in terms of the "unphysical" time variable t 

s{i) = ^{e+p + e-p), (45) 

I t^+p J I t^-p if o 3^ 1 

m = { , (46) 

lP + \ni ifp=l 

pP - 1 

m = ^5^^, (47) 

V{^) = Vo{l-^ni+^)f', (48) 

where a{p) = 2 + 1 /p and f3{p) = 2 — 1/ p. This class of solutions is only consistent 
with a flat FRW geometry. As we shall find in what follows, although it depends 
on the two parameters A and n, its behaviour is essentially determined by the n 
parameter. 

The cosmo logical density parameter for this class of solutions is ^2;^ = 1, while 
the deceleration parameter is given in terms of the "unphysical" time variable t 

,ii) = l^i'' + I^P'' + 'i^ .49) 

^^ ^ (n + p)H^P - g2t> + (n - p)2 ' ^ ^ 

where 

qi{n) = {n + p){l - n), (50) 

g2(n) = 2(/-n2), (51) 

g3(n) = (n-p)(l -n). (52) 

10 



The behaviour of the class of solutions (45)-(48) is illustrated in Figure 1 
where S{t) and V{iIj) are plotted for a representative set of values of the A,n 
parameters. 

Two points immediately stand out from eqs.(45) and (46). First, from eq.(46) 
we can see that for < n < 3 we have 

0<t<+oo =^ < t < +00, (53) 

while for n > 3 we have 

< t < +00 =^ - oo <t < +00. (54) 

Secondly, from eq.(45) we can see that all the solutions of this class corresponding 
to the range | < n < 3 are big-bang solutions, that is, start expanding from a 
Friedmann singularity (S* = 0), while those corresponding to the ranges < n < | 
and n > 3 are singularity- free. Indeed, for all solutions with | < n < 3 the edge 
S* = is always reached in a finite proper time into the past (S* ^ as t ^ 0) (see 
Fig. 1(c)). On the other hand, there is no Friedmann singularity for < n < | 
because these solutions always "bounce" at a finite proper time into the past 
given by (see Fig. 1(a)) 

_ ^ /p-n\'^p 1 /p-n\ ''^ 

neither for n > 3 because, in this case, an infinite proper time into the past is 
required to reach the edge S* = (see Fig. 1(e)). 

The asymptotic behaviour of the class of solutions (45)-(48) as t ^ is given 

by 

S{t)^{ " , (56) 



t™ 


if < n < 3 




1 {-tr 


if n > 3 




-1 + i^ot' 


if < n < 


3 


-i + M- 


-ty if n > 3 




\/(V')~ 


(1 + ^)^ 





^(t) ~ { , (57) 

(58) 
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where 

, s n — p 2p 

m^) = -, , s(n) = , (59) 

I — p 1 — p 

and 

^o = (2|l-p|)^ (60) 

the asymptotic behaviour of the n = 3 solution resuhs from the above equations 
in the hmit n — *> 3. Hence, we see that for an initial inflationary phase to 
occur we must have m > 1 4^ < n < |orn> 1. More precisely, the 
solutions corresponding to the range 1 < n < 3 exhibit an initial power-law 
inflationary phase of expansion (S* > 0, S* > and H < 0) which, in the limit 
n — > 3, yield to an initial phase of standard De Sitter exponential inflation; 
those corresponding to the range n > 3 display an initial phase of pole inflation 
(superinflationary expansion, S,S,H all positive); finally, those corresponding 
to the range < n < | exhibit an initial phase of contraction characterised by 
S* < 0, S* > and H > 0, which leads to a bounce. In this situation the criterion 
based on q to define inflation, does not hold, since S* < 0. 

It is worth noticing that while in the minimal coupling case superinflation 
only occurs for the k = +1 models [21], in the conformal coupling framework 
that can occur for other values of the curvature index, namely for k = 0. 

The late time behaviour (t -^ +oo) of the class of solutions (45)-(48) is given 

by 

S{t) ~ t\ (61) 

^(t)~l-^oot~^ (62) 

\/(^)~(l-^)^ (63) 



where 



and 



Kn) = —— , r n = -^, 64 

l+p I +p 



V'oo = (2(l+p))-^ (65) 

so, we see that for a late time power-law inflationary phase to occur we must have 
I > 1 <^ n > 1, the limit n -^ +oo corresponding to a late time standard De 
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Sitter exponential inflationary phase. While in the minimal coupling case power- 
law inflation is normally driven by exponential potentials V{iIj) oc exp(— A-?/'), 
with A constant > 0, we flnd from eqs.(61),(63) that in the conformal coupling 
case power-law inflation is driven by a polynomial potential. 

All this conclusions on the asymptotic behaviour of the class of solutions 
(45)-(48) are corroborated by a direct analysis of the asymptotic behaviour of 
the deceleration parameter q (cf. eq.(49)). Indeed, as t ^ g asymptotes to 

[n — pY n — p 

yielding q{t) < only when 0<n<|orn>l, while its late time behaviour 
(t -^ +oo) is given by 

g(t) ^ ^L^ = 1^^, (67) 

yielding q(t) < only when n > 1. Note that in these asymptotic inflationary 
behaviours the slow-rolling approximation (see Ref.[23]) is not valid. 

An exceptional case occurs when n = ^: on the one hand, the potential is flat 
V{iIj) = and on the other, we can write the corresponding solution explicitly in 
terms of the physical cosmic time t 



M = _l+.!^, (68) 

A 2 V 12 ' ^ ' 



As t — !► -|-oo the cosmic scale factor asymptotes to 

S{t) ~ f/\ (70) 

Thus, from a volumetric point of view, a flat FRW geometry fllled with a massless 
and non-interacting homogeneous conformally coupled scalar fleld exhibits as late 
time behaviour that of the standard flat FRW radiation-dominated solution. This 
should be compared with the analogous situation in the Einstein frame, where a 
homogeneous free scalar fleld propagating in a flat FRW background behaves as 

S{t) ~ f/\ (71) 

13 



To summarize all these results on the global and asymptotic behaviour of the 
class of solutions (45)-(48) we have collected them together in Table 1. 

4.2 Linear Expansion 

Next, consider the solution in the Einstein frame characterised by a linear 
expansion behaviour 

S{i) = At, (72) 

4!{i) = lniP, (73) 

V{^) = Voexp(--A, (74) 

where A is a strictly positive constant, 

//12 , r,\ 1/2 

P(k,A)=[^^j , (75) 

and Vo{p) = C V^? which is always possible for k = 0, +1 and can be made 
possible for k = —1 when subjected to the constraint A > 1. This class of 
coasting solutions may be viewed as a modern version of the Milne solution (which 
is characterised hj p = p = and k = —1 [17]), where the scalar field matter 
source allows it to be non-empty and with a generalization that any value of the 
curvature index k is allowed [16]. 

Mapping the class of solutions (72)-(74) in the Einstein frame into the physical 
frame with the help of eqs.(38)-(41), we obtain the following parametric class of 
solutions in terms of the "unphysical" time variable t 

Sii) = ^{i'^^ + i'~n, (76) 

m = { , (77) 

ip + ilnt ifp=l 

pP - 1 

m = ^5^r^, (78) 

l^(^) = \/o(l-^)"(l + ^)^ (79) 

14 



where a{p) = 2 + 1 /p and j3{p) = 2 — 1/p. This class of solutions also is consistent 
with all the FRW geometries, the A-amplitude being in the k = —1 case subjected 
to the same constraint as in the Einstein frame. 

The cosmological density parameter for this class of solutions is given in terms 
of the "unphysical" time variable t by 

'^^ K - k)i^P + {UJ2 - 2k)pP + (CU3 - k) ' ^ ^ 

where 

uJi{k,A) = A\l+py + k, (81) 

uj2{k,A) = 2A\l-p^) + 2k, (82) 

uj3{k,A) = A\l-pf + k; (83) 

and the deceleration parameter is given by 

The behaviour of the cosmic density parameter of the class of solutions (76)- 
(79), eq.(80), is illustrated in Figure 2 for closed and open FRW geometries. For 
each of these cases a representative set of values of the p parameter is plotted. 

Two points immediately stand out from eqs.(76) and (77). First, from eq.(77) 
we can see that for 0<p<l-^A'^> k/2 (this inequality is trivially satisfied 
by all the flat and open solutions of the class, being also satisfied by the closed 
solutions with A^ > |) we have 

0<i<+oo ^ 0<t<+oo (85) 

while for p > 1 ^^ A"^ < k/2 (this inequality is only satisfied by the closed 
solutions of the class with A^ < |) we have 

0<t<+oo ^ -oo < t < +00. (86) 

Secondly, from eq.(76) we find that all the solutions of this class corresponding 
to the range < p < 1 are big-bang solutions, while those corresponding to the 
range p > 1 are singularity- free. There is no Friedmann singularity for p > 1 
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because these solutions always "bounce" at a finite proper time into the past (see 
Fig. 3(a)), given by 



tbo 



p+i 



p-l\ ^P 



7) 



S{t) 



2(1 +p) \p+lj 2(1 -p) \p+lj ' 

neither for p = 1 because, in this case, the corresponding solution starts expanding 
from a non-zero value of the cosmic scale factor 6*00 (see Fig. 3(b)). 

The asymptotic behaviour of the class of solutions (76)-(79) as t ^ is given 

by 

t if p < 1 

Soo ifp=l , (88) 

(-t) ifj9>l 

' -l+ipot' iip<l 

-l+exp(4t) iip=l , (89) 

^ -1+M-ty ifp>l 

V{ij)r^{l+ij)f, (90) 



~ < 



m 



~ < 



where 



2p 



sip) - . 

I — p 

iJo = mi~p\y- 

while the late time behaviour (t -^ +00) is given by 

S{t) ~ t, 

^(t)~l-7/>oor^ 



where 



r(p) 



2p 



(91) 
(92) 

(93) 

(94) 
(95) 



(96) 
(97) 



1 — p 

V'oo = (2(l+p))-^ 

Hence, we see that the class of solutions (76)-(79) does not exhibit an asymptotic 
infiationary phase neither as t — ^ nor as t ^ +00 approaching instead, in both 
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cases and for general p, the S{t) ~ t behaviour of the coasting solutions. This 
conclusion is corroborated by a direct analysis of the asymptotic behaviour of the 
deceleration parameter since g ^ both as t ^ and as t ^ +00 (cf. eq.(84)). 
Although the class of solutions (76)- (79) is not asymptotically inflationary 
according to the negative deceleration parameter criterion, it has been argued that 
the coasting solutions ought to be considered inflationary [21]. In fact, if the linear 
expansion continues for long enough it will solve all the well known kinematical 
problems of the standard hot big-bang model. Furthermore, as shown in Ref. 
[22], these models allow the generation and evolution of density perturbations 
that may be responsible for the large scale structure that we observe in the 
Universe. However, unlike what is usually expected from inflationary models 
the present models can lead to ^2^ 7^ 1 today [21]. This interesting feature of the 
coasting models, which may provide a solution to the cosmic dark matter problem 
[23] , can be illustrated by the asymptotic behaviour of the cosmological density 
parameter corresponding to the class of solutions (76)-(79), since as t ^ we 
have (see Fig. 4) 



and as t ^ +00 we have 



"'<*' ~:^k-' + Anhw- <*' 



4.3 De Sitter Expansion From a Singularity 

Finally, let us consider the solution in the Einstein frame characterised by a 
"sinh" expansion from a Friedmann singularity (S* = 0) 

^(t) = Asinh(cjt), (100) 

m = In ('tanh^(^)') , (101) 

VW = Vo Q +/sinh=^(iv^)j , (102) 
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where A, oj are strictly positive constants, 

/I k V^ 

and Vq = C"^^^, which is always possible for /c = 0, +1 and can be made possible 
for k = —1 when subjected to the constraint Auj > 1. For large t, this class of 
solutions approaches a standard De Sitter inflationary phase, the cosmic scale 
factor growing exponentially with time in this limit. 

Mapping the class of solutions (100)-(102) in the Einstein frame into the 
physical frame with the help of eqs.(38)-(41), we obtain the following parametric 
class of solutions in terms of the "unphysical" time variable t 

S{i) = ^ siiih(u)i) ('tanh>'(y ) + coth''(y )") , (104) 

^ = i(tanh^(f)+coth^(f)). (105) 

^ tanh^^(f)-l_ 

^' tanh2P(f) + l ^ ' 

VW = V,{1- ^2)2 + V2 ((1 + tpni - i,f + (1 + i,f{l - ^r) , (107) 

where 

a{p) = 2 + - , /5(p) = 2--, (108) 

p p 

and 

V^{u;,A) = 2u'-j^ , V,{uj,A) = ^(^uj' + ^y (109) 

Here p is restricted to be a positive integer and A, uj are flxed in such a way as 
to satisfy the following constraint equation 

AV ^ 3-i-^, (110) 

Although the starting class of solutions in the Einstein frame eqs.(100)-(102) 
is consistent with all the FRW geometries, this class is only consistent with a 
closed FRW geometry. This reduction in the number of admissible FRW geo- 
metric backgrounds, is a consequence of the integration of eq.(105) only being 



1^ 



possible when p is a positive integer, which happens only when k = +1. How- 
ever, the integrals of integer powers of the "tanh" and "coth" functions are given 
inductively [24] so, we can not write the solution of eq.(105) in closed form for 
general p. To accomplish this we have to restrict ourselves to a particular value 
of p. 

In what follows we will focus, for simplicity, on the p = 1 solution. Taking 
p = 1 in eqs.(104)-(107), we get 

S{i) = A cosh{ui), (111) 

t(t) = -ln(sinh(cjt)), (112) 

CO 

^-\i) = -cosh{uji), (113) 

\/(V^) = \/o(l-^'), (114) 

where A^u^ = \ and Vq = 3/2A^. 

Now, it's interesting to note that eq.(112) can be solved for t, which enables 
one to eliminate the t parameter between eqs.(lll)-(113) and to write the solution 
explicitly in terms of the physical cosmic time t. When this is done, we obtain 

S{t) = A{1 + exp{2ujt)Y/^, (115) 

iP(t) = -(1 + exp(2cjt))-i/^ (116) 

with V{iIj) given by eq.(114). 

The cosmological density parameter for this solution is given in terms of the 
physical cosmic time t by 

Q^{t) = 1 + 2(exp(-2tut) + exp(-4cut)), (117) 

and the deceleration parameter is given by 

q{t) = - (^ + exp(-2cut)) . (118) 

For very large t {t ^ +oo), the cosmic scale factor asymptotes to 

5(t) ~exp(cut), (119) 
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while for t ^ — oo we have 

S{t) ~ A. (120) 

This is a non-singular inflationary solution that approaches asymptotically a stan- 
dard De Sitter exponential inflationary phase. Such late time inflationary phase 
results here, as in the minimal coupling case, from the presence of the constant 
term Vq in eq.(114), which is just an effective cosmological constant. On the 
other hand, the cosmological density parameter makes a smooth transition from 
infinity {Qd -^ -|-oo as t ^ — oo) to unity {Qd ^ 1 as t ^ +00), so the solution 
is asymptotically flat. 

It is worth noticing that the slow- rolling condition (see Ref. [23]) although 
asymptotically satisfied in the Einstein frame, fails to be true in the physical 
frame. The ratio A = lip/ipl/H is equal to unity instead of A ^ 1. So, to a 
slow-rolling inflationary phase in the Einstein frame do not always corresponds 
a slow-rolling inflationary phase in the original frame. In fact, the condition for 
this to happen requires ■?/'^ ^ 1 as shown in Ref. [8], which is violated in our case 

ir < 1). 

5 Summary and Conclusions 

In this work we have found a number of new, exact classes of conformal 
scalar field cosmologies. The solution-generating technique used here extends 
the technique originally introduced by Bekenstein for massless conformal scalar 
fields [14], and is based on the conformal equivalence between the non- minimal 
coupling model characterised by the action (1) and canonical Einstein's gravity 
plus a canonical scalar field. This technique has been used in the literature to 
map a number of familiar results from the Einstein frame, such as the reheating 
temperature and density spectrum, into the physical frame [25]. It has also been 
used to derive new exact cosmological solutions to some scalar-tensor theories 
[26]. Nevertheless, to our knowledge, the conformal transformation technique has 
never been used to derive self-interacting conformal scalar field cosmologies. 



20 



We have applied the extended solution-generating technique to the classes of 
solutions in the Einstein frame that were recently obtained by Ellis and Mad- 
sen [16]. For the classes of solutions obtained in the physical frame we have 
discussed the existence of singularities, the asymptotic behaviour and the pos- 
sibility of occurrence of inflationary phases. As a result, we have found that 
one class of solutions (eqs.(45)-(48)) tends asymptotically to a power-law infla- 
tionary behaviour S{t) ~ t^ with p > 1 yielding, in the limit p -^ +00, to a 
late time phase of standard De Sitter exponential inflation; and that another 
class (eqs.(76)-(79)) exhibits a late time approach to the S{t) ~ t behaviour of 
the coasting models, leading to any asymptotic value for the cosmological den- 
sity parameter, in particular to ^2^ 7^ 1, depending on the choice of the initial 
conditions. Both of these classes comprise bouncing solutions which avoid an 
initial singularity. It was also found that, a flat FRW geometry fllled with a 
massless and non-interacting homogeneous conformally coupled scalar fleld ap- 
proaches asymptotically the typical S{t) ~ t^/^ behaviour of the standard flat 
FRW radiation-dominated solution. Our results thus reveal that to simple be- 
haviours of the cosmic scale factor in the Einstein frame corresponds a wealth of 
behaviours in the physical frame, namely, we flnd several different phases, with 
in some cases, a smooth transition between ordinary and inflationary expansion. 
Note that to a slow-rolling inflationary phase in the Einstein frame do not always 
corresponds a slow-rolling inflationary phase in the physical frame. This indicates 
that some caution is required in transposing to the physical frame results which 
are derived in the Einstein frame, in particular when regarding inflation. 

This work also serves the purpose of establishing the classes of potentials in 
the physical frame which drive the various behaviours of the cosmic scale factor 
considered in the Einstein frame. This allows a comparison to be made with the 
corresponding potentials in the Einstein frame. For instance, while in the minimal 
coupling case a power-law inflationary phase is normally driven by exponential 
potentials V{iIj) oc exp(— A-?/'), with A constant > 0, we have found that in the 
conformal coupling case it is driven by a polynomial potential(cf. eq.(63)). 

Although we have focused on some classes of solutions in the Einstein frame. 
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the technique used here can also be utihsed to generate other conformal scalar field 
cosmologies by starting from other minimal scalar field cosmologies. Furthermore, 
the extented solution-generating technique can also be applied to other geometric 
backgrounds. In a forthcoming paper we will apply it to a number of anisotropic 
scalar field cosmologies. 
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Table Captions 

Table 1: Summary of the global and asymptotic volumetric behaviour of 
the solutions of the power-law expansion class (eqs.(45)-(48)). 



Figure Captions 

Figure l:The figure shows the scale factor plotted against the time t and 
the potential plotted against ip for the power-law expansion class (eqs.(45)-(48)): 
(a)-(b) n = I and A = 2 (^ p = y^), (c)-(d) n = | and A = 2 {^ p = y^), 
(e)-(f) n = 12 and A = 2 (^ p = 2). The graphs (a),(c) and (e) correspond to 
the scale factor. The graphs (b),(d) and (f) correspond to the potential. 

Figure 2:The figure shows the time evolution of the cosmological density 
parameter Q^ for the linear expansion class (eqs.(76)-(79)): (a) k = +1 and 
p = I (^ A = V3), (b) k = +1 andp=l {^ A = -^), (c) A; = +1 and p = 2 
i^A= ^), (d) k = -landp=^i^A = |). 



Figure 3:The figure shows the scale factor plotted against the time t for 
the solutions of the linear expansion class (eqs.(76)-(79)) corresponding to the 
parameter values (a) p = 2 (^ A = Jjj) and (h) p = I (^ A = 

Figure 4:The figure shows the asymptotic values taken by the cosmological 
density parameter ^2^ as t ^ (Qq) and as t ^ +oo (Qoo) plotted against the 
A parameter, for the solutions of the linear expansion class (eqs.(76)-(79)). The 
graphs (a) and (b) correspond to the closed {k = +1) model. In (a) the graph 
diverges to infinity as A ^ w|, asymptotes to 4 as A ^ and to 1 as A ^ +oo. 
The graphs (c) and (d) correspond to the open {k = —1) model. 
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Table 1 





Global 
Behaviour 


Asymptotic Behaviour 


i^O 


i^ +00 


< n < 1/3 


Bouncing Models 


Contraction 


Ordinary Expansion 


1/3 < n < 1 


"Standard" Big-Bang Models 


Ordinary Expansion 


Ordinary Expansion 


n= 1 


"Standard" Big-Bang Model 


Linear Expansion 


Linear Expansion 


1 < n < 3 


Inflationary Big-Bang Models 


Power-Law Inflation 


Power-Law Inflation 


n = 3 


Non-Singular Inflationary Model 


De Sitter Expansion 


Power-Law Inflation 


n > 3 


Non-Singular Inflationary Models 


Pole Inflation 


Power-Law Inflation 
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